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INTRODUCTION 


In  attempting  telemetry  of  data  from  test  bodies  re-entering 
the  earth's  atmosphere,  the  problem  of  penetrating  with  an  electro¬ 
magnetic  wave  the  sheath  of  ionized  gases  (plasma)  which  envelopes 
these  bodies  at  hypersonic  velocities  is  encountered.  A  knowledge  of 
reflection  and  transmission  coefficients  for  the  sheath  is  of  value  in 
attempts  at  solution  of  this  problem. 

This  paper  presents  a  method  for  computing  reflection  and 
transmission  coefficients  for  a  plane -parallel  inhomogeneous  isotropic 
layer  of  plasma  when  the  inhomogeneity  is  a  function  only  of  distance 
along  a  normal  to  the  surface  of  the  layer.  A  plane  wave  at  normal  in¬ 
cidence  is  assumed. 

It  is  recognized  that  as  a  model  for  the  plasma  sheath  described 
above,  the  idealization  treated  is  deficient  in  several  important  respects. 
It  neglects  the  curvature  of  the  layer  as  well  as  the  fact  that  the  incident 
wave  is  not  plane.  The  effects  of  induction  fields  near  the  antenna  and 
the  problem  of  antenna  breakdown  are  similarly  ignored. 

Despite  these  objections,  it  is  believed  that  solutions  of  the 
simplified  problem  may  be  taken  as  order  of  magnitude  estimates  of  the 
transparencies  of  sheaths  having  the  prescribed  distributions  of  electrical 
properties.  Moreover,  the  results  should  be  of  interest  for  comparison 
with  those  of  more  comprehensive  calculations. 
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This  and  related  problems  have  been  treated  by  several  authors, 

but  most  of  these  neglect  the  effect  of  spatial  variation  of  collision  fre- 
7 

quency  which  Kritz  has  shown  to  be  significant  (in  certain  cases)  in 
determining  the  reflection  and  transmission  characteristics  of  the  plasma. 

Exact  analytic  solutions  in  closed  form  are  possible  only  for  the 
simpler  distributions  of  permittivity  and  conductivity,  and  more  com¬ 
plicated  distributions  have  been  treated  using  various  approximation 
1,  2,  6,  8 

methods  .  Existing  solutions  of  these  types  for  a  few  simple 

3 

plasma  geometries  are  discussed  and  tabulated  by  Graf  and  Bachynski  . 

These  solutions,  however,  have  rather  limited  applicability  and  may  prove 

cumbersome  in  calculations  of  reflection  and  transmission  coefficients. 

A  direct  and  expedient  method  for  computing  these  coefficients  for  a  wide 

variety  of  distributions  of  both  permittivity  and  conductivity  which  takes 

advantage  of  available  automatic  computation  facilities  is  needed. 

These  requirements  are  probably  best  met  by  a  method  involving 

the  direct  numerical  integration  of  Maxwell's  equations  within  the 

plasma.  This  approach  is  not  new,  and  several  workers  in  the  field 

have  discussed  the  relative  merits  of  different  forms  for  the  equations 

7 

and  of  various  integration  processes.  Kritz  presents  a  method,  and 

„  .  10  .  t  „ 

Zivanovic  uses  a  numerical  technique  to  compute  a  set  of  matrix 

elements  descriptive  of  the  plasma  layer  after  the  manner  of  four  terminal 

networks  and  transmission  lines  in  network  theory.  Klein^  and  Budden* 


2 


treat  the  problem  of  joining  approximate  analytic  solutions  with  numerical 

solutions  at  boundaries  separating  the  regions  of  validity  of  the  two. 

Each  of  these  approaches  offers  certain  advantages. 

The  computer  program  used  provides  flexibility  in  the  choice  of 

distributions  of  both  conductivity  and  permittivity  and  the  modified 

Runge-Kutta  integration  process  employed  allows  for  control  of  accuracy 

in  the  solution.  The  program  is  described  in  greater  detail  in  a  report 

5 

by  Kavanaugh  and  Scarborough  . 
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ANALYSIS 


Maxwell's  equations  for  a  stationary  medium  containing  no  free 
charges  are: 


V  •  D(t) 

V  •  B(t) 

V  x  E(t) 

V  X  H(t) 


0 

0 


J  + 


3D(t) 

0t 


(U 


For  an  inhomogeneous  isotropic  plasma,  the  constitutive  equa¬ 
tions  may  be  written 

B(t)  =  pH(t) 


D(t)  =  e  E(t) 


(2) 


J  =  <rE(t) 

with  p  ,  c.  and  <r  scalar  functions  of  position. 

In  a  plasma  it  is  observed  that  the  effective  permeability 

is  very  nearly  equal  to  the  free -space  value  pQ.  It  has  been 

5,  6,  8, 

shown  in  several  works  that,  if  only  those  fields  whose  time 

-  iu)t 

dependence  may  be  expressed  as  e  are  considered,  the  conductivity 
effective  in  the  plasma  is  given  by  the  relation 


vnet 


<r  = 


T(V2  +  CJ2) 


(3) 
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This  expression  is  arrived  at  by  considering  the  motion  of  an 

electron  under  the  influence  of  an  electromagnetic  wave  and  subject  to 

— ►  — ► 

damping  by  collisions  with  heavier  particles.  The  v  x  B(t)  term  in  the 
Lorentz  force  is  found  to  be  negligible. 

The  same  analysis  leads  to  the  following  expression  for  the 
permittivity  effective  in  the  plasma: 


e  =  e 


1  - 


ne 


me  (v2  4-  to2) 
o 


(4) 


ne 


In  terms  of  the  plasma  frequency  wp  the  expressions  (3) 

and  (4)  are 

... .  2 


o  (v*  +  <j2) 


(5) 


and 


e  =  e 


1  - 


v2  +  w2 


(6) 


-*■  -*•  -iu)t  .  “*  -  icjt 

Writing  E(t)  *  E  e  and  H(t)  m  H  e  and  using  the 


constitutive  equations  (2),  Maxwell's  equations  (1)  become 


V  •  e  E  =  0 


V-  H  =  0 

V  X  E  =  iwu  H 

o 

V  X  H  =  (<r  -  iwe)  E  . 


(7) 
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Confining  attention  to  the  particular  geometry  of  interest,  if  the 
spatial  variations  of  the  field  quantities  and  e  and  ( r  depend  on  a  single 
coordinate  z  normal  to  the  surface  of  the  layer,  the  harmonic  fields 
are  transverse  and  the  equations  (7)  reduce  to 


3  c  E 
3z 


=  0 


-  3H 
n  •  - —  =  0 
3  z 

-  3E  - 

n  X  - —  =  iu ix  H 
3  z  o 


(8) 


3H 
n  X  - — 
3  z 


uu  ,  v  zr 

—  =  (or  -  iwe)  E 


in  which  n  is  a  unit  vector  along  the  z  axis.  Expanding  the  first  of  these 
yields  the  relation 

-  -•  d  e  0E 

n  '  E  TT  +  £  n  •  Tz  =  0  •  (9) 

The  first  term  of  (9)  represents  a  coupling  between  the  electric  field 
vector  and  the  gradient  of  the  inhomogeneity  in  permittivity  which  is  here 
zero  since  n  .  E  =  0  for  a  transverse  field  and  normal  incidence. 

Hence,  the  first  of  equations  (8)  becomes 


3E 

3z 


=  0  . 
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A  further  simplification  of  the  equations  is  possible,  since  only 
the  magnitudes  of  the  field  need  be  considered,  their  directions  in  space 
being  constant.  The  third  and  fourth  of  equations  (8)  then  reduce  to 


dE 

dz 


i  oj  u  H 

O 


dH 

dz 


i  (e  co  +  i  a)  E  > 


(10) 


the  first  two  equations  yielding  no  additional  information  concerning  the 
fields. 


Equations  (10)  have  been  solved  exactly  for  only  a  few  simple  dis¬ 
tributions,  t  (z)  and  <r(z)  and  the  general  case  must  be  treated  by  various 
approximation  methods  or  by  numerical  techniques.  Since  in  the 
present  instance,  numerical  values  for  reflection  and  transmission  co¬ 
efficients  are  of  greater  interest  than  a  general  solution  to  (10)  and 
facilities  for  high-speed  automatic  computation  are  available,  the  latter 
alternative  is  favored. 

7 

Following  a  procedure  suggested  by  the  work  of  Kritz  and 
Budden*,  a  transmitted  wave  E^  of  a  particular  amplitude  (unity)  and 
phase  (kQz  -  k0zQ),  is  assumed  in  the  free  space  region  immediately 
’’outside"  (z  >  zQ)  the  plasma  layer.  (See  Figure  1.)  With  (z  )  =  1 


and  H  1  (zq)  = 


V 


as  initial  values,  the  four  simultaneous  equations 


represented  by  (10)  are  integrated  numerically  over  the  region  0  <  z<  zQ, 
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HOMOGENEOUS 

DIELECTRIC 


e  p  u  ,  cr  =  0 
m  o 


tt1  xr  1  ikz 
E  =  E  e 

o 


H1  =  H  1  eikz 
o 


R  _  R  -ikz 
E  =  E  e 
o 


HR  =  H  R  e*ikZ 
o 


ANALYTIC 


PLASMA 


e  (z).  [iQ,  o-  (z) 


dE  u 

—  =  1U)LL  H 

dz  o 


=  i  (t  uj  +  i  o")  E 
dz 


NUMERICAL 


Figure  1.  Schematic  Representation  Showing  the 
Type  of  Solution  Valid  in  Each  Region 


8 


the  integration  proceeding  backwards  along  the  z  axis.  At  the  boundary 
z  =  0,  it  is  required  that  the  E  and  H  fields  be  continuous,  i.  e. ,  that 


E1  (0)  +  ER  (0)  =  E  (0) 
H1  (0)  +  HR  (0)  =  H  (0) 


(ID 


But 


E1  (0)  =  E  1  , 

o 

er  (0)  =  e  r 

O 


(0)  =  /J^o1* 


(12) 


and 


HR  (0) 


-21  R  ‘ 
1*  0 


Substituting  from  (12)  into  (11),  the  following  pair  of  simultaneous 
equations  is  obtained: 


E  1  +  E  R  =  E  (0) 

o  O 


E  X-E  R 
o  o 


(13) 


—  H  (0) 
m 


Solving  for  E0^  and  E0^  gives 

and 


-j 

E  (0)  +  /  H  (0) 

m 


(14) 


E 


R 


o 


E  (0)  - 


9 


The  reflection  coefficient,  defined  by 


R  = 


(15) 


i  E 


I,2 


becomes,  in  terms  of  Er,  E^  ,  Hr,  Hj  defined  by  E(o)  2  Er  +  i  Ei  and 
H  (o)  «  Hr  +  iHp 


R  = 


(JT  E  -  n/jx  H  )  +  ('s/e  E.  -  n/ul  H.)J 

m  r  or _ m  x  ro  1 

(n/T~  E  +  <s/jr  H  )Z  +  (s[T~  E.  +  JT  H.)2 
m  r  o  r  mi  *o  1 


(16) 


Similarly,  the  transmission  coefficient  is 


T  = 


c  E 
o  1  o 


4  \fe  ~ 
m  o 


C  I  r*  1  I2 

m  |Eq  | 


(\ ft~  E  +  n/JT  H  )2  +  (\ir~  E.  +  n/JT  H.)2 
m  r  or  mi  o  1 


(17) 


Requiring  conservation  of  energy  gives  for  the  absorption  coefficient 


A  =  1  -  (R  +  T) 


(18) 


U 


EXAMPLE 


The  following  study  will  serve  as  an  example  of  the  type  problem 
to  which  the  method  is  applicable: 

It  is  required  to  find  transmission  and  reflection  coefficients  for 
the  plasma  sheath  surrounding  certain  bodies  during  re-entry  at  an  alti¬ 
tude  of  80  km  and  at  speeds  of  Mach  18  and  Mach  30.  Various 
positions  on  the  body  surface  are  to  be  considered  in  order  to  determine 
the  most  favorable  location  for  a  transmitting  slot  antenna,  and  the  effect 
on  transmission  of  varying  the  dielectric  constant  of  the  "window"  cover¬ 
ing  the  slot  is  to  be  investigated.  A  transmitting  frequency  of  240  Mc/s 
is  assumed. 

Electrical  properties  of  the  plasma  were  computed  from  chemical 
equilibrium  flow  field  values  of  pressure  and  enthalpy.  Two  points  on  the 
body  surface  (at  2  Rn  and  4  Rn,  as  measured  along  the  axis  of  symmetry, 
where  Rn  is  the  nose  radius  of  the  body)  were  selected  as  representative, 
(Figure  2).  Profiles  of  relative  permittivity  and  conductivity  along  lin^s 
normal  to  the  body  surface  at  these  points  were  taken  as  the  distribuv  *ns 
of  electrical  properties  within  a  plane -parallel  plasma  slab  and  along  a 
normal  to  its  surface.  A  more  complete  description  of  the  above  proce¬ 
dure  is  to  be  found  in  the  appendix,  together  with  curves  representing  the 
distributions  prevailing  at  the  various  speeds  and  positions. 

The  distributions  were  approximated  by  the  following  analytical 
expressions: 


11 


Geometry  of  Re-entry  Body 


Mach  18 


2R 


1  -  0.  4  e 


-108  z 


0.  04  e 


-38. 4  z 


4R 


''l  -  0.09  e-'56'8  Z 


0.  01 35  e 


-28.  3  z 


Mach  30 


JR 


1R , 


■8.  V  ,  /  *  0.85 


104  5  .  -  I  0  I 


-  =  *■?  :  0.534 


0  <  z  <  0 .  0 95 

z  >  0.  095 

0  <  z  <  0.  095 

z  >  0.  095 

0  <  z  <  0.  095 

z  >  0.  095 

0  <  z  <  0.  095 

z  >  0.  095 


U  <:  z  <0.  0681 
z  >  0. 0681 

0  <  z  <  0.  0980 

z  >  0. 0980 

0  <  0. 0976 

z  >  0. 976 

0  <  7.  <  0.  0954 

z  >  0. 0954 


0.  1  m 


t*  -  0.  1508X10 10  radian /sec 


RESULTS 


The  reflection  and  transmission  coefficients  computed  using 
these  distributions  and  for  several  values  of  em  are  tabulated  below: 


Ma.ch  30 

2  Rn  4  Rn 


•  m/£  o 

R 

T 

R 

T 

-  3 

-  3 

1 

0.  998 

0. 529  x 

10 

0.998 

0.  478  x  10 

2 

0.  998 

0.742  x 

io’3 

0.997 

0.669  x  10‘3 

3 

0.  997 

0.  903  x 

10~3 

0.997 

0. 810  x  10"5 

-2 

-  3 

4 

0.  997 

0.  104  x 

10 

0.  997 

0.925  x  10 

-2 

-2 

5 

0.997 

0. 115  x 

10 

0.996 

0.  102  x  10 

Mach  18 

2  Rn 

4  Rn 

Til*  O 

R 

T 

R 

T 

1 

0. 989  x 

10'4 

0.995 

0.  175  x 

10“4 

0.998 

2 

0. 290  X 

10-1 

0.996 

0. 293  x 

io'1 

0.968 

.  1 

-  1 

3 

0.711  x 

10 

0.  924 

0.716  x 

10 

0.926 

4 

o.  no  x 

10° 

0.  885 

0. Ill  x 

10° 

0.  887 

5 

0.  145  x 

10° 

0.  851 

0.  146  x 

10° 

0.852 

Reflection  and  Transmission  Coefficients  for  the  Plasma  Sheath  Sur¬ 
rounding  a  Body  Re-entering  the  Atmosphere  at  Mach  18  and  Mach  30, 
at  Two  Positions  on  the  Surface  and  for  Various  Dielectric  Constants 
of  the  Antenna  Window. 
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CONCLUSIONS 


Several  conclusions  of  a  qualitative  nature  may  be  drawn  from 
these  results.  At  Mach  18,  transmission  is  almost  complete  with  very 
little  reflection  occurring,  whereas,  at  Mach  30,  the  situation  is  re¬ 
versed.  It  will  be  noted  that  in  neither  case  is  absorption  of  energy 
responsible  for  appreciable  loss  in  transmission,  since  the  absorption 
coefficient  A  *  1  -  (R  +-  T)  is  in  every  case  very  small.  The  effect  is 

due  almost  entirely  to  increased  reflection. 

Also  of  interest  is  the  effect  on  transmission  of  the  value  of 
permittivity  of  the  ’’window’*.  At  Mach  18,  an  increase  in  this  per¬ 
mittivity  results  in  a  decrease  in  transmitted  power,  whereas  the  opposite 
is  true  at  Mach  30.  It  is  suggested  that  further  investigation  of  this 
effect  is  needed  to  determine  whether  it  is  significant.  The  results 
further  indicate  that  at  Mach  30,  absorption  is  more  pronounced  at  the 
4Rn  position  than  at  2Rn,  as  evidenced  by  the  smaller  values  for  both 
R  and  T  at  this  position.  The  reason  for  this  is  not  yet  apparent,  but 
the  difference  in  wake  thickness,  together  with  the  fact  that  c  has  large 
negative  values  at  this  speed,  are  probably  responsible. 

Summarizing,  the  results  indicate  that  for  the  assumed  trans¬ 
mission  frequency,  failure  to  penetrate  the  Mach  18  sheath  may  not  be 
attributed  to  reflection  or  absorption  of  the  electromagnetic  energy  in  the 
wave  by  the  plasma  through  the  mechanisms  considered  here.  At  Mach  30, 


15 


however,  these  effects  almost  certainly  will  preclude  transmission 


at  this  frequency  and  at  any  position  on  the  body  surface. 
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APPENDIX 


Values  of  pressure  and  enthalpy  in  the  flow  field  surrounding  a 

particular  body  were  computed  assuming  chemical  equilibrium  for  speeds 

of  Mach  18  and  Mach  30  at  an  altitude  of  80  km.  ^ 

From  these  values,  the  corresponding  temperatures  and  mass 

density  ratios  were  determined  from  a  Mollier  diagram  of  properties  of 
(a) 

equilibrium  air.  Electron  densities  were  then  read  from  a  second 
chart  giving  electron  density  as  a  function  of  mass  density  ratio  for 
various  temperatures.  ^  Values  for  the  collision  frequency  v  were 
determined  from  the  formula 


v 


/i  t/v  v  pressure  in  atmospheres 

l.  10  x  10‘5  x  - - - : — rzr- - 

temperature  in  K 


(c) 


Using  the  values  of  n  and  v  thus  determined,  <r  and  t  were  com¬ 
puted  using  equations  (3)  and  (4).  Contours  of  constant  c  /e0  are  plotted 
in  Figures  3  and  5,  and  contours  of  constant  <r/e0w  are  plotted  in  Figures 
4  and  6,  and  approximate  location  of  the  shock  wave  is  also  shown. 

Plot 8  of  t  /  t0  and  °7*0u>  as  functions  of  distance  along  the 
normals  are  given  in  Figures  7  through  10.  The  Mach  18  curves  are 
approximated  by  exponential  functions;  the  Mach  30  curves  by  linear  func¬ 
tions.  These  analytic  approximations  appear  in  the  main  body  of  the  report. 
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